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Abstract 



A family of one-dimensional quantum systems described by certain weakly 
non-Hermitian Hamiltonians with real spectra is studied. The presence of 
a microscopic spatial defect (i.e., of a fundamental length-scale 9 > 0) is 
mimicked by the replacement of the real line of coordinates 1 by a non-tree 
graph. By construction, this graph only differs from R on a small interval 
of size O(0). An unexpected return of the reality of the whole spectrum in 
anomalous strong-coupling regime (i.e., the emergence of an isolated island 
of stability of the system) is then revealed, proved and tentatively attributed 
to the topological nontriviality of the graph. 



1 Introduction 

1.1 Schrodinger equations on linear lattices 

For any concrete Hamiltonian with real spectrum and with the property H ^ 
valid in an auxiliary, unphysical Hilbert space one must introduce 
another, "standard" physical Hilbert space in a way explained, say, 

in Ref. [I]. Most often this is achieved via the physics-determining metric 
operator = 0^ ^ I which must be compatible with the "cryptohermiticty" 
condition 

H = q-^ H ]q.- H \ (1) 
guaranteeing the necessary Hermiticity of H in the correct Hilbert space 

Perturbation-expansion constructions of the required ad hoc Hilbert-space 
metrics = Q(H) are difficult (cf. their characteristic sample in Ref. [3]). 
In what follows the simpler approach will be employed in which one approx- 
imates the (real) interval of coordinates by a discrete lattice composed, say, 
of N = 2K grid points 

, £-K+2 , • • • , £-2 , £-1 , £o , £l , £ 2 , • • • , £,K-1 , £k ■ (2) 

This converts the solution of Eq. ([1]) to the straightforward application of 
linear algebra. 

Without additional interactions the free motion confined to lattice fl2]) 
with K < 00 may be controlled by the Runge-Kutta recipe leading to the 
2K— dimensional matrix Schrodinger equation 

(3) 

A systematic amendment of precision may be mediated by an ad hoc increase 
of K. For more details the readers may consult our recent paper jl] on bound 
states (to be abbreviated as BI in what follows). Lattice (J2]) with finite 
K < 00 has been interpreted there as an approximate representation of the 
straight real line. Equally well one can turn attention to the unbounded 
equidistant grid with K = 00 (cf., e.g., our papers |5J on scattering). 

In paper BI we tried to enhance the phenomenological appeal of the model 
and introduced a specific VT— symmetric interaction in it. In a way inspired 
by Ref. [S] the purely kinetic Hamiltonian of Eq. has been replaced there 
by the manifestly non-Hermitian tridiagonal and IK— dimensional matrix 
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we extracted the respective spectra of energies E = En (u),n = 1, 2, 
N = 2K. Empirically we revealed that these spectra remain strictly real 
inside a K— independent weak-coupling interval of v G (—1, 1). 

The key problem addressed in paper BI was the construction of the com- 
plete family of metrics or, equivalently, of all the eligible inner products 
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1p*(£m) @m,n0(£n) 



(4) 



between elements ip and <fi of the "standard" Hilbert space H^ s \ We sim- 
plified this construction by working with the non-complex, real Hamiltonian 
matrices and with their real left and right eigenvectors. Then, the resulting 
metrics (which can be, in general, complex-valued) were also real. 

Even then, the selection of the appropriate (i.e., symmetric and positive- 
definite) matrix = 0(H) in Eq. @ was not unique []]. Naturally, the 
ambiguity of O represents one of the main weaknesses of the theory. Among 
its suppressions available in the literature one may mention not only the most 
common (often called Dirac's) unit-matrix choice of Q( Dirac ) = / but also the 
increasingly popular specification of metric 0^ cv ' = CV were V is parity 
while symbol C denotes a charge (review paper [6] may be recommended as a 
source of more details). In this context the main physical message delivered 
by our paper BI may be read as the conjecture of the classification of the 
eligible physical metrics in terms of a certain microscopic fundamental- 
length scale 9 G M + . Within such a scheme one can treat 0^ Dirac ' and Q^ cv > 
as extreme special cases characterized by 6 = and 9 = oo, respectively. 



1.2 Nonlinear lattices and graphs simulating the finite- 
size microscopic defects 

One of the first applications of the freedom of the choice of any finite 9 G 
(0, oo) was described in our paper BII [7]. We tried to move there beyond 
the one-dimensional quantum systems and, tentatively, we replaced the real 
line of coordinates M by the q— pointed star-shaped graph G^ q > embedded, 
presumably, in a more- dimensional space M. d . 



Fortunately, the discretization method proved applicable at all the not 
too large integers q > 2. At the same time, the scaling-invariance of the 
corresponding tree-shaped quantum graphs G^ created a subtle inconsis- 
tency in the new theory. Indeed, on a pragmatic phenomenological level the 
scaling invariance of graphs G^ contradicts the existence of a finite (though, 
presumably, very small) "smearing length" 9 < oo removing the ambiguity 
of the physical metric G. 

From another perspective, the presence of a non-vanishing "smearing 
length" 9 implies the admissibility of topologically nontrivial anomalies in 
the graphs themselves. In what follows we intend to develop this idea in 
more detail. In one of the simplest realizations of such a selfconsistently 
non-local quantum-graph project we shall replace the real line of paper BI 
by a suitable short-range graph-supported modification of the real line. 



Figure 1: The simplest topologically nontrivial localized modification of the 
real interval of coordinates. 

For the sake of simplicity we shall only contemplate the simplest non-tree 
graph sampled in Figure HJ with the admissible size of the defect G — R 
restricted just to the order of magnitude of 9. 



2 Non- Hermit ian interactions near vertices 
2.1 Single- loop model 

There exist many topologically nontrivial (i.e., non-tree, loops-containing) 
generalizations of the trivial real-line graph (i.e., of R = G^ in our present 
notation) which might be discretized along the lines discussed in papers BI 
and BII. Here, in a way guided by Fig. [TJwe shall select the first nontrivial, 
non-tree discrete-graph lattice in the form 
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(5) 

possessing four vertices X-k, x\ and xk and four wedges. For the sake of 
simplicity of our present constructive considerations just the external wedges 
will be assumed of variable length, K = 1,2,.... 



In the next step following the current practice [8] we shall start from the 
use of the most common discrete Laplacean A on this lattice, with 
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where we choose w = 2, and with 
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where we choose u — 3. 

In the final step the natural generalization of the model will be obtained 
when we append elementary Hermiticity-violating nearest-neighbor interac- 
tion terms to the two outmost vertices X-k and xk and, independently, also 
to the remaining two inner vertices x_i and x±. For illustration purposes 
we may take graph fl5]) and indicate the position of the interaction terms by 
symbols 4k (representing the elements proportional to a real coupling g), £ 
(representing the coupling h) and 9? (representing the coupling z). At any 
K this will define our present three-parametric family of toy Hamiltonians 
H = H^ K \g, h; z). At an illustrative choice of K = 3 this means that we 
may pick up the one-loop discrete graph 
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(6) 

as leading to the eight-dimensional sparse-matrix Hamiltonian 

H^(g,h;z)= (7) 
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2.2 Factorizable secular equation 

From the purely computational point of view a remarkable property of our 
K = 3 model is that its secular equation gets factorized. Indeed, once we put 
H = H^i^i + S, 7 — 5; z) we reveal that the energies may be identified with 
the roots of one of the following pair of the two algebraic (and, in principle, 
exactly solvable) polynomial secular equations 

E 4 - 9 E 3 + P± E 2 + Q± E + R± = (8) 



with the respective coefficients 

P+ = P+(z, 7) = z 2 + 24 + 4 7 2 , Q + = Q + {z, 7) = -5 z 2 - 19 - 16 7 2 , 

i?+ = i2+(z,7) = 2z 2 + 47 2 2 2 + 12 7 2 + 2 (9) 
(which do not depend on 5) and 

p_ = p_( z ,6) = 28 + z 2 + 45 2 , Q_ = Q_(z,5) = -35 - 5 z 2 - 16 5 2 , 

R- = R„(z,5) = U + Qz 2 + 125 2 + A5 2 z 2 (10) 

(which do not depend on 7). This means that at each "outer" coupling z the 
spectrum is composed of the two one-parametric quadruplets of levels which 
might be expressed in closed form, in principle at least. 
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Figure 2: Typical spectral pattern: Eight K = 3 energies as functions of z 
at fixed g = h = 0. 

A numerical illustration of the parametric dependence of the spectrum is 
provided by Fig. |2j We see there that the whole octuplet of the z— dependent 
energies remains real if and only if we stay in the weak-coupling regime with 
\z\ < 1. The first merger occurs at \z\ = 1 involving the second and the third 
root of the + — subscripted secular polynomial flS}. In the strong-coupling 
domain, i.e., at \z\ > 1 these two energies form a complex-conjugate pair 
with the growing common function |Imi? + |. 

The separation of the "physical" weak-coupling regime from its "unphys- 
ical" strong-coupling complement is "generic" within a fairly broad class of 
models. This has been demonstrated in our purely numerical study [9] where 
the "generic" spectrum has been found composed of a part called "robust" 
(i.e., real at all the coupling strengths) and a complement called "fragile" 
(sampled by the two ellipses in Fig. [2]). The same or similar pattern has also 
been detected via the exactly solvable models [TU] (offering also the explicit 
forms of the imaginary parts of the energies) as well as via their simplest, 
semi- numerically tractable generalizations [TT] . 

In such a setting our present introduction of topologically nontrivial quan- 
tum graphs may be characterized as opening new perspectives and moving 
beyond the classification offered by Ref. [9]. 

After the change of notation with g = 7 + 5 and h = 7 — 5 we may keep 
5 = and let the second parameter grow, 7 = g = h > 0. Empirically we 
reveal (and, subsequently, easily prove by elementary means) that the overall 
pattern represented by Fig. [2] gets only inessentially deformed in the whole 



interval of 7 G (0,1). For illustration, Fig. [3] displays the shape (i.e., the 
smoothly deformed form) of the spectrum at g = h — 0.98 (i.e., near the 
weak-coupling physical horizon). 




Figure 3: The spectrum sampled by Fig. [2] at 7 = gets almost doubly 
degenerate at 7 = 0.98. 



3 Anomalous domains of stability 

The inspection of coefficients ([9]) and ( TTTJT) in secular polynomials (JSj) reveals 
that the contribution of the growth of the positive quantity 7 2 or S 2 shares 
the direction with the contribution of z 2 . Hence, these pairs of parameters 
play a partially interchangeable dynamical role. No particularly surprising 
changes of the overall form of the spectrum can be expected at small 7, 
therefore. Moreover, from the purely qualitative point of view one of the 
three parameters can be considered redundant so that we shall simplify the 
discussion by setting 5 = in what follows. 

3.1 Numerical experiments with Hamiltonian iif (^(7, 7; z) 
in the strong-coupling regime 

In the specific K = 3 model with vanishing 5 = g — h = and subcritical 
positive 7 < 1 our numerical experiments indicated that all the spectrum 
may become doubly degenerate in the limit 7 — > 1. Due to the relative 
simplicity of the model this expectation can be rigorously confirmed by the 
explicit evaluation of the secular eigenvalue condition at 7 = 1, 

(E-2) 2 [E 3 ~7E 2 + (14 + z 2 ) E -7 - 3 z 2 ] 2 = 0. (11) 

One may now ask what happens when one moves beyond this "natural" 
boundary of the weak-coupling cryptohermiticity domain, i.e., in the lan- 
guage of mathematics, beyond the boundary of the well-established domain 
of the guaranteed reality of the spectrum, i.e., from the point of view of 
quantum physics, beyond the boundary of the domain of the safely stable 
and safely unitary time-evolution of the system. 

A deeper change of the pattern can be expected to emerge in the upper 
vicinity of the critical parameters, i.e., say, at 7 > 1. We performed a number 
of numerical experiments in such a new domain, with a typical result obtained 
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Figure 4: The atypical spectral pattern in the strong-coupling regime with 
7 = 1.035. 

at 7 = 1.035 and sampled in Fig. |H In full accord with expectations one 
reveals, first of all, that the spectrum ceases to be realm the half-strong/half- 
weak-coupling regime with \z\ < 1. Indeed, the numerical spectrum becomes 
unphysical there since a pair of the energies becomes complex (this "missing 
pair" is absent in the picture of course). In parallel one can notice that 
during the growth of 7 > 1 the two lowest energy levels still stay very close 
to each other forming an almost degenerate doublet while the constant root 
E = 2 ceases to be degenerate. 

The comparison of Figs. [3] and H] at \z\ > 1 (i.e., in the genuine strong- 
coupling regime) reveals that while the maximal (i.e., the eighth) energy 
remains safely real, the seventh and sixth level merge and complexify beyond 
certain fairly large pair of "exceptional-point" values of z — ±\z( EP > 
(numerically one finds that 21^^(7) ~ ±3 at 7 = 1.035). 
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Figure 5: The fine-tuned reality of the spectrum in the strong-coupling 
regime: At 7 = 1.035, the zeros of the secular determinant stay all real 
for z G (1, 1 + a) with a ~ 0.022. 

For these reasons (easily supported also by an elementary algebra) there 
still exists a comparatively smaller domain of strong couplings where the 
situation remains ambiguous. A magnified detail of Fig. H] displayed in Fig. |5] 
reveals the presence of an extremely anomalous behavior of the numerical 
energy spectrum at 7 > 1 and \z\ > 1. Indeed, while the two largest energies 
(lying out of the frame of this picture) stay safely real we discover that all 
the remaining sextuplet of the energies also remains real inside the interval 
ofz G (1,1 + a) w (1.000,1.022). 

Such an anomalous return of the complete stability and observability of 
the system in the strong-coupling regime looks surprising and unexpected. 



Certainly, the confirmation of its existence must be performed by non-numerical 
means. Fortunately, such a proof is rendered possible by the factorizability 
(jHD of the underlying secular equation. 



3.2 The domain of reality of the first, 7— independent 
quadruplet of energies 

The study of the spectrum of Hamiltonian if ^(7, 7; z) (where 6 = 0) gets 
perceivably facilitated by the 7— independence of the minus-subscripted co- 
efficients (ITU]) in one of our quartic polynomial secular equations (jHD- 

The quadruplet of the related (here, minus-superscripted) energies is easy 
to find because one of them is constant, E$ = 2 while the other three roots 
will only vary with z in a way dictated by their definition 

-7 -3z 2 + UE + E z 2 -7 E 2 + E 3 = . 

At \z\ = 1 this equation may be easily factorized so that the triplet of ener- 
gies Ej\z), j = 1,2,3 will coincide with the explicitly known roots of the 
polynomial 

Q(E) = -10 + 15 E -7 E 2 + E 3 = (E -2) ^-L±^ 

Once we consider positive z 2 — 1 = A, the graphical representation of our 
secular equation in the form of the sum Q(E) + (E — 3) A = indicates that 
with the (unlimited) growth of A G (0, 00) the rightmost root E$ \z) will 
stay real. It will merely decrease from its maximum E l f \±l) ~ 3.618 to 
its minimum E { 3 \±oo) = 3. Hence, we may restrict our discussion to the 
specification of the domain of the survival of the reality of the two smaller 
roots e[~2(z). 

There certainly exists an interval of A G (0, \ m ax) in which the latter roots 
remain real while moving towards each other. This motion starts at their 
respective initial values e[ (±1) ~ 1.382 and E% (±1) = 2 and terminates 
at the common degenerate value of 

E[-\±y/l + X max ) = E ( -\±^l + \ max ) = y. 

The latter point is the position of an extreme (in fact, of the maximum) of our 
secular polynomial so that its value must satisfy the conditions Q'(y) + A = 
and Q"(y) < 0. The latter one is safely satisfied since we surely have y < 7/3. 

The former condition looks tractable as quadratic equation with the two 
roots y of which we need the smaller one, 



7 - a/4 - 3 A 
V = 3^ • 

Nevertheless, for our present purposes we should rather treat this condition 
differently, viz., as the linear definition of A = X(y). The insertion of this 
quantity in secular equation finally yields the single cubic equation for y, 

-35 + 42y- 16y 2 + 2y 3 = 0. 



This equation gives the unique real value y ~ 1.702843492 of the maxi- 
mal energy which agrees with the prediction displayed in Fig. The ex- 
act analytic formula also exists for the unique maximal admissible value of 
^max ~ 0.14078102. Its knowledge provides the first rigorous necessary con- 
dition of validity of our numerically supported stability-pattern predictions. 



3.3 The domain of reality of the second, 7— dependent 
part of the spectrum 

A deeper inspection of Fig. Vindicates that the plus-subscripted secular equa- 
tion OH]) + (EJ) should be expected more restrictive. Still, its analysis may 
proceed along similar lines, using positive A = z 2 — 1 and fi = 7 2 — 1 in our 
second secular equation 

E 4 - 9 E 3 + (29 + A + 4 //) E 2 + 

+ (-40 - 16// - 5 A) E + 16 + 2 A + 4 (1 + //) (1 + A) + 12 fj, = . 
At the boundary A = fi = the left-hand-side polynomial gets factorized, 

Q{E) = E A - 9 E 3 + 29 E 2 - 40 E + 20 = (E 2 - 5 E + 5) (E - 2) 2 . 

In the strong-coupling regime with A > and fi > the graphical interpre- 
tation of the plus-subscripted secular equation may again be based on the 
split of the secular polynomial into the easily factorizable polynomial Q(E) 
and correction term 

R(E) = E 2 \ + 4 E 2 \i - 16 Efi - 5 E\ + 6 A + 16 fj. + 4 fx A 

exhibiting the three-component form 

R(E) =4\li + \(E-2) (E -3) +Afi(E - 2) 2 . 

Obviously, polynomial Q{E) vanishes at its four roots E^ ~ 1.382, E^ = 
E^ = 2 and E^ ~ 3.618 so that the spectrum is real at A = fi = 0. 
Now, our task is to prove that this quadruplet of energies remains real in a 
non-empty domain T> of nontrivial A > and [i > 0. 
After the reparametrization of our variables 

E= x+^^ /i = i6/i, A = 4A 

we get the simpler form of the secular equation, 

S(X, jl, x) = (x 2 + fi- 5) (x + l) 2 + A (x 2 + fi - 1) = (12) 
which may be further rewritten in the non-polynomial form 

A 4 

(13) 



(x + 1) 2 x 2 + fi-l 

comparing the left-hand-side expression Flhs[{x + I) 2 , A] (which is solely 
controlled by parameter A and exhibits the reflection symmetry with respect 



to the off-central point x = — 1) with the right-hand-side symmetric function 
F RHS {x 2 , p) of x (which is independent of A). 

Obviously, the latter equation cannot have any real roots for p G (5, oo) 
because in this interval we have F LH s[(x + l) 2 , A] > 1 while F RHS (x 2 , p) < 1. 
With the decrease of p < 5 there emerges, on the purely geometric grounds, 
the first degenerate and positive pair of real roots £3,4 (A, p) which lies on the 
right branch of the spike Flhs[(x + I) 2 , A]. With the further decrease of p 
(and at any A > 0) the rightmost root 2:4 (A, p) can only move to the right 
so that its reality remains granted. Similarly, its neighbor £3 (A, fx) moves to 
the left and stays also real and bigger than — 1. 

This means that the reality of the spectrum is controlled by the reality 
of the pair of the smaller roots £1,2 (A, P) which have to lie on the left branch 
of the spike F LHS [(x + l) 2 , A]. The existence of the latter two intersections 
with the right-hand-side curve Frhs{% 2 , P) requires that p < 4. As long as 
the domain of stability T> must lie in the strip of not too large jx G (0, 4) 
with bounded A G (0, X max (p)), let us now return to the polynomial version 
of Eq. (fT2"j) . paying attention just to its relevant x < — 1 part. 

In this interval the reality of the two leftmost roots £1.2 (A, p) of poly- 
nomial S(X,p,x) becomes more easily interpreted after we separate this 
polynomial in two components, viz., S(X,p,x) = Q(p,x) + XR(p,x) where 
Q(p, x) = (x 2 + p — 5) (x + l) 2 . The reason is that R(p, x) = x 2 + p — 1 is 
safely positive in all the interval of x < — 1. On this background the reality 
of the roots xi^{X,p) may be deduced from the facts that 

• the auxiliary A = curve Q(p,x) stays negative (and has a minimum) 
in non-empty interval of x G (—6, —1) where b = —Xi(0,p) = a/5 — jx; 

• secular polynomial S(X,p,x) stays negative (and has a minimum at 
some point y) inside a smaller interval of x G {—b + A 1; —1 — A 2 ); 

• both the shifts Ai^ vanish at A = and grow with growing A < 

Xmax (A) • 

The auxiliary quantity y = y{p) and the function A = X max (p) which de- 
termines the boundary of the domain T> are both defined by the pair of 
polynomial equations 

S(X max (p),p,y) = 0, S'(X ma x{P),P,y) = 0. 

The easy elimination of X max (p) yields the final implicit exact definition of 
the unknown auxiliary parameter y = y(p), 

+ 2 py 2 - 2y 2 + 4y + p 2 - 6 p + 5 = . 

We may re-read this relation as the explicit definition of the doublet of the 
eligible inverse functions 

P = P±(V) = -V 2 + 3 ± 2 y/-y2+l-y. (14) 

They remain real for y G (— c, — 1) with c = (1 + v / 5)/2 ~ 1.618. Using 
Taylor series it is easy to prove that also the smaller function p-(y) remains 



positive in the right vicinity of y = —c. The inspection of Fig. [6] indicates, 
moreover, that the lower bound of fi-(y) > —0.4 is still not too negative so 

that with g = h = \Jl + yu._ (t/) / 16 > 0.987 we may keep speaking about the 
strong coupling regime. 



Figure 6: The pair of auxiliary functions fi±(y). 




Figure 7: The pair of bracketing functions X ma x{y)- 

By the backward insertion we obtain the explicit form of the doublet of 
the functions of our main interest, viz., of 



(y + 1) (y 2 + y - 2 T 2 ^-y 2 + I - y) 



^maxlKv)] = Kuxlfiiv)] = 7— j • ( 15 ) 

Eqs. ( fl4|) and (115p form our final result. They provide the parametric defini- 
tion of the boundary of the domain T> where the anomalous, strong-coupling 
stability of the quantum model with Hamiltonian H^( , y, 7; z) takes place. 
Thus, any parameter y E (— c, — 1) bounded by c = (1 + Vo)/2 ~ 1.618 
specifies the coupling constants g = h = yT+ fi±(y)/lQ via Eq. ( !T4l) . The 
spectrum of energies will be then real for all the couplings z £ (1, z max ) re- 



stricted by the upper bound z max = y 1 + A^ a:r (?/)/4 which is determined by 
formula (fl5l) . 



4 Discussion 



As long as our secular polynomial of the eighth degree gets factorized into 
two polynomials (IE|) of the fourth degree, one could, in principle, employ 



the discriminant techniques for quartic polynomials and try to extract the 
reality conditions for the spectrum. In fact, our present analytic approach 
to the description of the stability islands in parameter space represents just 
a formally equivalent, purely algebraic implementation of such an idea. 

In doing so we heavily profitted from additional simplifications offered by 
our particular example in a way recommended and used in a similar study of a 
certain specific VT— symmetric chain model family in Ref. [12] • In this man- 
ner, the unbelievably clumsy explicit discriminants (or, alternatively, their 
merely perturbative versions) were replaced by a perceivably more efficient 
equivalent specification of the exceptional points. 

This being said it is necessary to add that the main weak point of such 
an algebraic approach to the problem of horizons of the physical domains lies 
in its insensitivity to perturbations of various kinds. In particular, also our 
present determination of the anomalous behavior of the energies as sampled 
by Figure |5] may prove sensitive to perturbations. The simplest possible 
form of such a perturbation may be mimicked by a trivial scalar shift of the 
secular determinant making it either slightly positive or slightly negative. 
Such a type of a purely numerical experiment reveals the occurrence of the 
two alternative scenarios where the complexification of the energies caused 
by the decrease of z > 1 involves either the mere two central energy levels 
(cf. Fig. [8]) or the two non-central neighboring energy-level-pairs (cf. Fig. [9]), 
respectively. 
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Figure 8: A change of the spectrum of Figure [5] under a perturbation. 




Figure 9: A change of the spectrum of Figure [5] under another perturbation. 



A more systematic classification of the similar alternatives would already 
move us beyond the scope of our present paper. Interested readers could 



find a complementary study of this problem in Ref. [13] where a purely com- 
binatorial full classification of all of the possible confluences of the energy 
levels (i.e., of a complete menu of eligible quantum catastrophes) has been 
found, in full generality, due to the extreme simplicity of the underlying toy 
model. In an opposite extreme, Ref. [14] could be consulted for a realistic 
physical example where the authors were able to identify certain helical tur- 
bulence function and its derivative as effective parameters responsible for the 
different kinds of the unfoldings of eigenvalues in a complicated magnetohy- 
drodynamic system. 

For a deeper formal insight in the theory the Kato's monograph [15] could 
be consulted. Supplementary comments on the practical aspects of the ex- 
plicit determination of the boundaries dT> (carrying, usually, the nickname of 
"exceptional points" [12]) may be sought not only in magnetohydrodynam- 
ics [T7] but also in the context of classical mechanics [IB] or experimental 
optics [T9] . 

Recently, the very real correspondence between the exceptional points 
and stability of quantum systems attracted attention to the study of con- 
crete models. In many of them the necessary assumption of Hermiticity of 
observables takes place "in disguise" . For illustration of details one can rec- 
ollect, e.g., the studies of the interacting boson models in nuclear physics 
[1] as well as of perturbation expansions in field theory [20]. Many new ap- 
plications as reviewed, e.g., by Bender [6] range from the supersymmetric 
systems [21 J to integrable models [22] and from the first-quantized relativis- 
ts systems [23] to the model-building in quantum cosmology [24] . Important 
innovations also occurred in the description of the time-dependent quantum 
systems [25] or of the theory of scattering [5j [26] . 

In many of these applications an important role is played by the proofs 
of the reality of the spectrum. Often, it can only be mediated by perturba- 
tion theory [27] so that the difference between our present "normal" , weak- 
coupling real spectrum and its 'anomalous" strong-coupling version might 
be perceived as purely conventional. Similarly, the role of the "smallness 
of size" of our present anomalous domain of parameters guaranteeing the 
stability of the system may be also declared just supportive. This being 
said, we still believe that the ultimate source of the emergence of an iso- 
lated island of stability in our present model should be sought in the implicit 
continuous-limit connection of this model with the topologically anomalous, 
scaling-non-invariant graph of Fig. [TJ 

5 Summary 

In a toy model a loss of locality on microscopic distances has been simulated 
by the brute-force replacement of the real line of coordinates by a topologi- 
cally nontrivial loop-containing graph. We also recalled our recent papers BI 
and BII and complemented such a "kinematical" form of the loss of locality 
by its parallel, dynamically generated form in which the interaction ceases 
to be naively Hermitian. 



Both our toy-model kinematic and dynamical nonlocalities were assumed 
restricted to a small-size spatial domain. On this background we performed 
certain numerical experiments which revealed certain spectral anomalies in- 
volving, i.a., the puzzling disappearance of the instabilities in the strong- 
coupling dynamical regime. The existence of this phenomenon has been con- 
firmed by its thorough non-numerical description and analytic explanation. 

The resulting anomalous parametric-dependence of the spectra certainly 
reflects the presence of the hypothetical microscopic deviation of the space in 
which we live from its standard topologically trivial picture. We believe that 
the similar quantum phenomena indicating possible deeper connections be- 
tween kinematical topology and dynamical analysis will deserve an intensive 
study in the future. Our present encouraging technical message is that the 
use of the quantum graphs and of nontrivial discrete lattices seems to offer 
one of the feasible mathematical ways towards the similar phenomenological 
analyses. 
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